Testing the reliability of a velocity definition in dispersive medium 
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We introduce a method to test if a given velocity definition corresponds to an actual physical flow 
in a dispersive medium. We utilize the equivalence of the pulse dynamics in the real-cj and real-fc 
Fourier expansion approaches as a test tool. To demonstrate our method, we take the definition 
introduced by Peatross et al. [Phys. Rev. Lett. 84, 2370 (2000)] and calculate the velocity in 
two different ways. We calculate i) the mean arrival time between two positions in space, using the 
real-o; Fourier expansion for the fields and ii) the mean spatial displacement between two points in 
time, using the Fourier expansion in real-fc space. We compare the velocities calculated in the two 
approaches. If the velocity definition truly corresponds to an actual flow, the two velocities must 
be the same. However, we show that the two velocities differ significantly (3%) in the region of 
superluminal propagation even for the successful definition of Peatross et al. 

PACS numbers: 42.25.Bs, 41.20.Jb 



I. INTRODUCTION 

Studies concerning the propagation of light in disper- 
sive media dates back as far as BriUouin and Sommerfeld 
Nevertheless, new studies on the concept of 

pulse velocity were stimulated by the famous experiment 
[isl [iGj , where light seems to propagate over the speed of 
light in vacuum (superluminal). This effect takes place 
due to superluminal group velocities near the absorption 
resonance in dye solutions. 

Beside absorptive dielectrics there exist metamateri- 
als where index exhibit some unusual behaviour. These 
are constructed either using coherent population trap- 
ping e.g. Electromagnetically induced transparency 
(EIT) and index enhancement schemes, or spatial mod- 
ulation as in left handed materials [l^. Beyond scien- 
tific curiosity, applications such as memory storage/read 
out using slow light motion carried out the physics of 
metamaterials to an important status in the information 
technology. 

In such materials, dielectric response may switch from 
absorption to gain in MHz intervals for the optical fre- 
quencies. As well, index of refraction may display steep 
and negative derivatives. In those situations, it is con- 
fusing even to figure out the reorganization of the pulse 
beside facing the problem of superluminality. Unfortu- 
nately, when the original pulse is severely modified there 
is no direct way to test the validity of the propagation 
velocity. The experiments 0-@ measure either the peak 
of the pulse or mean absorption time. However, motion 
of the pulse peak or center may not correspond to a travel 
velocity since the shape of the pulse is distorted by mu- 
tual act of gain/absorption. For this reason, we adopt a 
method to test the reliability for a given description of 
velocity. We utilize the equivalence of the pulse dynam- 
ics within the real-w and real-fc Fourier expansions as the 
test tool in complex dielectrics. 

When the dielectric function is complex, there 
are two alternative mathematical approaches that 



are used in dealing with pulse dynamics. One 
can analyse the system equivalently using the real- 
ty {E{x,t) = /+^da;£>i(a;)e^('=(")^-"*)) and real-fc 

{E{x,t) = /+^dfci:»2(fc)e'(''''-"('='*)) Fourier domain. 
E{x,t) is the Electric Field that is used to calculate the 
velocity. In the real-w or real-fc approaches, decay ap- 
pears spatially or temporarily in the Fourier integrand, 
respectively, since fc or w is complex. If one is interested 
in the penetration depth, Fourier expansion is carried 
out over the rcal-cj space. On the other hand, if one 
is interested in the temporal lifetime of the pulse in the 
material, then real-fc expansion is used. In example, deal- 
ing with photonic crystals composed of complex dielectric 
materials, one constructs the Master equation using real- 
ty; (real-fc) Fourier space for calculating the penetration 
(duration) of light into (in) the crystal |19l - [2l| . For a 
given velocity definition, pulse speed can be calculated 
within both approaches. 

Among different velocity definitions 0, [13, existing 
in the literature one of the most successful is the one that 
is introduced by Peatross et al. Q- In Ref. [1], propa- 
gation is described with Poynting vector average of the 
temporal position, i.e. {t)x = J dttS{x,t)/ J dtS{x,t). 
The velocity vi — {x — Xo)/{{t)x — {t)xo) is introduced 
in this regard. On the other hand, considering the same 
definition Q, a second velocity W2 = {{x) t — {x) tg) / {t — to) 
can be adopted similarly using the mean spatial position 
of the pulse, i.e. {x)t = J dx x S(x,t) / J dxS{x,t). Since 
t or X average is dealt within the calculation of vi or 
f2, respectively, it is standard to work in the conjugate 
Fourier space where w or fc is chosen as real. If the def- 
inition correctly addresses a physical flow, then the two 
velocities must be identical or at least must be very close. 

It is shown [ll| that the observed consistency of the 
definition Q with the experimental results follows 
from the equivalence of the detector time (mean time for 
detector absorption) to the arrival time deduced from 
this description Q. Accordingly, here we choose to test 
the validity of this definition as the example. 
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In order to compare the two results for the given veloc- 
ity definition, we perform the following. We first calcu- 
late the mean arrival time At = {t)x-{t)xo, from position 
xq to X, for a distance Axi = x — xq. This is handled 
in the real-w approach. Second, we calculate the cor- 
responding mean propagation distance Ax2 ~ {x)t-(x)t„, 
from time to to t, in the real-fc Fourier expansion. For the 
purpose of comparison, we chose t — to equal to At which 
is the value determined in the real-w expansion. After- 
wards, we compare the two distances, Axi and Ax2, for 
the same At. Thereby, we compare the two velocities 
vi = Axi/At and V2 = Ax2/At. 

Paper is organized as follows. In Sec. II, we estab- 
lish a connection between the two expansion coefficients 
Di{lu) and D2{k) using the boundary conditions. In Sec. 
Ill, we calculate the velocity definition of Peatross et al. 
in two different ways, by expanding the fields both in 
the real-w and rcal-fc Fourier space. We show that there 
occurs differences in the calculated values of the two ve- 
locities, especially in the superluminal region. Sec. IV 
includes our conclusions. 




FIG. 1: (color online) Real-oj (Ci: red thin cj/ = line) 
and real-fc (C2: tliick contour) integration paths for index 
n{(jj) plotted for dielectric constant with Up = O.luo and 
7 — 0.12tJo- C2 contour is deduced from the nonlinear index 
equation ck — n(a;)cj for the real k values. The two lines Li 
and 1/2 are the branch cuts of the index [Sj. Both branch cuts 
are below the real-fc integration path. Length of the branch 
cuts are exaggerated only for the visual purposes. 



II. RELATING THE FOURIER COEFFICIENTS 

Di{lo) and Diik) 

In this section, we establish an analytical connection 
between the Fourier components of the real-w expansion 
{Di{uj)) and the real-A: expansion {D2{k)). In order to in- 
dicate that a variable is fixed to real, we use a bar accent 
over that variable. For instance, dCo (dk) corresponds to 
an integration over the real lu (k) space. 

We consider a dielectric function in the Lorentzian 
form El in 



e(w) = 1 



(1) 



where ujp (wo) is the plasma (atomic transition) frequency 
and 7 is the damping rate. One can calculate refractive 
index by n{uj) = (e)^^^- 

The calculation of complex k values for given real u) 
value is straightforward using ck = cjn((D). However, 
the calculation of complex w requires the solution of the 
nonlinear equation ck = ujn{uj) for given k real value. In 
Fig. [T] we plot both integration paths in the complex lo 
plane, corresponding to real-oj (Ci line) and real-fc (C2 
contour). We also indicate the branch cuts (^1,2) of index 
n{uj) I2J], which are below the rcal-fc integration path C2. 

In order to relate Fourier coefficients Di{ui) and D2{k), 
we consider the following simple reflection/transmission 
boundary problem (see Fig. [5]) . A Gaussian wave packet 
(of Fourier coefficient Ai{u!)) travelling towards right in 
vacuum {n = 1) is incident on the absorbing dielectric 
slab of index n{uj), see Fig. ^jji.. It results in a reffectcd 
wave packet of Fourier coefficient Bi(uj) and a transmit- 



ted wave (into the slab) 



+00 



Ei{x,t)= rfwDi(c^)e*('=(")^-"*) (2) 



with Fourier coefficient Di{uj). In Fig.[2}D, the same prob- 
lem is considered in the real-fc Fourier domain; with in- 
cident wave of coefficient A2{k), reflected wave of coeffi- 
cient ,62 (fc), and a transmitted wave 



E2ix,t) 



-j-00 



dkD2ik)e'^''''-'^^''^'^ (3) 



of Fourier coefficient D2{k). 

Since both lo and k are real on the left hand side, 
they arc simply related by A2{k) = cAi{uj) and B2{k) = 
cB2{uj). In the real-w approach, using the boundary con- 
dition (BC) at a; = 0, one obtains 



i +n(w) 



Ai(w). 



(4) 



In this paper, for the sake of simplicity, we consider a 
Gaussian profile for Ai{uj) which does not have any pole. 

Additionally, the two solutions for the transmitted 
wave, Ei{x,t) and E2{x,t), must match at the bound- 
ary x — for all times. That is 



+00 



dujDi{uj)e 



dkD2{k)e-"^^'^^\ (5) 



where w, k stand for real variables and uj{k) ^ cJi?(fc) + 
iujj{k) is complex function of the real variable k deter- 
mined from the nonlinear index equation ck = ujn{uj). 
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a) real-u! n = 1 


n = nji + in J 


x = 


b) real-k n = I 


n = njj + in/ 



X = 



FIG. 2: The same reflection/transmission problem considered 
(a) in the real-cj and (b) in the real-fc Fourier spaces. Incident 
light penetrates from vacuum on the left hand side (LHS) to 
an absorbing dielectric slab of complex index n{u}) on right 
hand side, (a) Incident wave of Fourier coefficients Ai(cj) 
results a reflected wave of coefficient Bi{uj) and a transmit- 
ted wave of coefficients Di{ijj). (b) Incident wave of Fourier 
coefficient A2{k) results a reflected wave of coefficient B2(k) 
and a transmitted wave of coefficient D2{k). Since on the 
LHS index is real and unity, it is simply A2{k) = cAi{uj) and 

B2{k)=cB2{Lu). 



Integrating both sides with j'^'^ dte^^^^'^ in Eq. 
one obtains 

-oc 

dkD2{k) / rfie^WC^ )-"('=))*. (6) 



+ 00 



Time integrations in Eq. ^ are in the form of Dirac- 
delta function (^(z)) that is generalized to complex ar- 
gument im 



/ + 00 
dte-'^^ = 2<(z), 
-GO 



(7) 



with the sampling property [26l | 

+cxD+'iIin{a} 

dzF(z)^(z - a) = F(a), (8) 

— oo+jlm{a} 

where a and z are complex. Thus, time integrals in Eq. 
© result ^ (w(fc') — d))) and ^ (a;(fc') — a;(fc))) , respec- 
tively. Now, Eq. becomes 



du]Di{uj)S,{uj{k') - Co) 



4-00 



dkD2{k)^{Lo{k')-L0{k)) , (9) 



where (D is real and uj{k) is complex. Note that dk corre- 
sponds to the 6*2 contour in Fig. [5] and all possible values 
of both Loik') and Loik) rely on the Co curve. 

Since and ^(aj(fc')— w) [231 have no pole or 

branch cut between Ci line and w/ = Im{w(fc')} line 
for any choice of k' (Fig. [5]), integration path in the LHS 
of Eq. ^ can be as well written as 



+cjo+Im{aj(/e')} 



dujDi{uj)£, {uj{k') - w)) = Di{uj{k)). (10) 

On the other hand, using dk = ^lo^ transformation 
RHS of Eq. © can be written as 



dkD2ik)^ {uj{k') - uj{k))) 



dk 



C2 



dLo—D2{k{<^))i{uj{k')-uj), (11) 
duj 



where C2 is the correspondence of the real-fc integra- 
tion path in the complex uj plane, depicted in Fig. [2] 
dk/duj = n{uj) + ujti'{uj) has no pole or branch cut in 
between the two contours Ci and C2. In addition, we 
assume that D2{k{Lj)) does not have pole or branch cut 
in between Ci and C2. Afterwards, we confirm this as- 
sumption. Therefore, we shift the C2 integration in Eq. 
([TT|) to the Ci contour and obtain 

+ °° rlh _ dk - - 

duj—D2{k{uj))^ {iu{k') ~cu) = —{k')D2{k'). 



(12) 

In obtaining the last term of Eq. (|12l) , in a similar way to 
Eq. pO)) . we further assume that D2{k{uj)) has no pole 
or branch cut between C'l line and w/ = Im{cj(fc')} line 
for any choice of k' . 

Finally, Eq. ([5]) transforms to the simple relation 



D2Ck) 



du) 
dk 



(fc)i?l(L.(fc)), 



(13) 



where presence of k indicates that all quantities are eval- 
uated at complex variable uj(k) corresponding to real 
fc. Eq. ([T^ confirms that D2{k) indeed does not have 
pole or branch cut in between the contours Ci, C2 and 
LOj = Im{aj(fc)} line. 

Relation (|13p is quite simple and straightforward. On 
the other hand, however, result is cumbersome when 
Di(u)) or n{Lu) has a branchcut/pole in between the 
curves Ci and C2 (see Fig. [1]). 

In the following section, we use Eq. and relate the 
real-fc integrand to the real-w one. 



III. COMPARISON OF VELOCITIES IN THE 
REAL-o; AND THE REAL-fc APPROACHES 

In this section, we test the velocity definition intro- 
duced by Peatross et al. 0. We derive the two expres- 
sions for the velocity, using the real-w and real-fc Fourier 
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FIG. 3: (a) Comparison of the two velocities (in units of c) 
deduced from definition of Peatross. et al. The velocities 
vi = Ax/ [{t)Ax — {t)o) and ?;2 = ((x)Ai — (a;)o)/Af are calcu- 
lated using the same definition, but performing real-cj (solid 
line) and real-fc (dashed line) Fourier expansions for the fields, 
respectively. For a consistent definition, the two results must 
be identical. However, 3% discrepancy between vi and V2 is 
observed in the superluminal region. Thus, this description is 
not so reliable in the superluminal regime. We use a Gaussian 
pulse of carrier frequency ujc and temporal width luot = 20. 
Propagation distance is Axi = 150c/a;o. (b) Real (nu) and 
imaginary (n/) parts of the index of refraction. Parameters 
for the index are the same with Fig. [1] 



expansions. Afterwards, we calculate the velocities vi 
and V2 by relating the coefficients I?i(w) and 1^2 (fc) (see 
Eg.s (|2l3p ) using Eg. (|T3| . In the superluminal region, we 
observe discrepancies in the amount of 3% (see Fig. [3^) ■ 

In Sec. II we derive the relation ()13[) imagining a 
boundary between vacuum and dielectric along the lines 
of experiments and the generality. In this section, 
however, we consider a more particular dispersive 
dielectric occupying the whole space without any bound- 
aries. This is done so, not to deal with the tail of the 
pulse relying out of the dielectric. Relations and 
remain valid. Because, the condition of matching the 
two solutions at the origin Ei{0,t) = E2{0,t) does not 
require a physical boundary. At any random position this 
condition must already be satisfied. 

We consider a Gaussian wave packet U{0,t) = 
e~* cos{ujct) imposed at the origin. This leads 



[23 to the Fourier coefficient Di{lu) = e"("""<= 



74 



where Wc is the carrier frequency of the 



pulse. I?2(fc) is determined from Di{uj) using Eq. ([Tc 

First, we calculate the arrival time {At = {t)Axi — 
{t)o) of the mean pulse center from position to Axi. 
Time averages are directly calculated within the real-w 
expansion of the fields similar to Ref . Q , using Fourier 
coefficient Di{u;). Second, we use the same arrival time 



At (which is calculated using real-w domain) in the real- 
k approach and evaluate the displacement of the average 
pulse position, i.e. Ax2 ~ (a;) At ~ (2^)01 from time to 
At. Finally, since At is common in both approaches, we 
compare the velocities wi = Axi/At and V2 = Ax2/At 
plotted in Fig. [3K- 



A. Real-oj 

Average time-position of the pulse after propagating a 
distance Ax (starting from x = 0) 



it) 



Ax 



J dttS{Ax,t) 
J dtS{Ax,t) ' 



(14) 



can be directly calculated by Fourier expansion in real-w 
space using Q 

dttSiAx,t) ^ Ax / dCj—e-'^'"^''\Di{Ld)fn*{u!) 
J-00 dw 



-\-oo 



duj 



and 



dtS{Ax, t) 



+ 00 



du)e-^'"^'=\Di{Co)\^n*{u)). (16) 



Here ki{u)) is the imaginary part of the wave- vector cor- 
responding to the real w value. The calculated values of 
the velocity vi ~ Axi/ {{t) ^^xi — for different carrier 
frequencies Wc are plotted in Fig. [3li with solid line. We 
choose a temporal width of lo^t ~ 20 and propagation 
distance of Aa;i = ISOc/wq. 



B. Real-A; 

Above, using the real-w approach, we determine the 
arrival time At in between the two positions and Axi. 
In the real-Zc approach, we use the calculated value of At 
as the input. We determine the distance that mean pulse 
center travels from time to At, i.e. Ax2 = {x)At — {x)o- 
The average pulse position at time At 



(x) 



At 



J dx X S{x, At) 
J dxS{x, At) 



(17) 



can be directly calculated by carrying the Fourier expan- 
sion over the rcal-fc coefficients using the expressions 

j dxxS{x,At) = At J^^ dk^e^'^'^' \D2{k)\'^ n*{k) 
+*/+^dfce2-^^*^i?2*(fc)n*(fc)(18) 



and 



dxS{x, At) 



+00 



dke 



2u}iAt 



D2{k)\n*{k), (19) 
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where Fourier components decay in time with imaginary 
part (w/) of the complex frequency ui{k) during the prop- 
agation. Complex uj{k) values are determined from the 
nonlinear index equation ck = un{u)) for real k, and uj 
is always negative for the absorbing dielectric Rel- 
evance of Fourier coefficient D2{k) to Di{oj) is given in 
Eq. (HH). 

The average displacement Ax2 that is calculated in the 
real-fc approach is compared with the one for the real-w 
approach Axi. We note that At is common to both ap- 
proaches. The calculated velocity U2 = [(x) At — (x) o) / 
is plotted in Fig. [5^ with dotted line for different carrier 
frequencies. The two results, vi and V2, differ signifi- 
cantly (3%) in the superluminal propagation regime. 

On the other hand, similar calculations using the real 
part of the conventional group velocity [21|, as vi = 
Kc{dLj/dk} and V2 = 1/Re {d/c/dw}, results in 16% dis- 
crepancy in the superluminal propagation region. 



IV. SUMMARY AND CONCLUSIONS 

The velocity introduced by keeping track of the pulse 
peak or the pulse center does not always corresponds 
to the velocity of the energy/signal transfer. When the 
pulse shape is modified during the propagation, it is con- 
fusing even conceptually to define the arrival time of the 
original signal. 

Here, we introduce a method to check if a given veloc- 



ity definition is reliable regarding its correspondence to a 
real physical flow. We calculate the velocity introduced 
by Peatross et al. [1] in two different ways. First, we cal- 
culate the mean arrival time At of the pulse between two 
points in space. We perform this calculation using real- 
cj Fourier expansion of the fields. Second, we calculate 
the mean displacement of the pulse between two points 
in time, and At. This calculation is carried out with 
rcal-fc Fourier expansion. Finally, since At is common in 
the both approaches, wc compare the two velocities. 

We observe that the velocity definition of Peatross et 
al., relying on the Poynting vector average of the pulse, 
results in 3% discrepancy in the superluminal propaga- 
tion region, see Fig.|3K. Thus, one questions if this veloc- 
ity truly corresponds to a physical flow in the superlumi- 
nal region. On the other hand, definition of Peatross et 
al. is more successful compared to the conventional defi- 
nition of group velocity, where discrepancy comes out to 
be 16% in the superluminal region. 

Since the arrival time introduced by Peatross et al. 
also corresponds to the detector time [ll|, we reach the 
additional conclusion that the arrival time measurements 
0-01 do not address a proper velocity for the flow. It is 
still an open problem to find a reliable velocity descrip- 
tion, consistent with equivalence of the two approaches. 
On the other hand, our method is also possible to ad- 
dress the physics of elementary particles when there exist 
sources standing for absorption/gain [28| . 



L. Brillouin, Wave Propagation and Group Velocity (Aca- 
demic Press, New York, 1960). 

A.I. Talukder, T. Haruta, and M. Tomita, Phys. Rev. 
Lett. 94, 223901 (2005). 

T. Kolimoto, Y. Fukui, S. Furue, K. Nakayama, and Y. 
Fukuda, Phys. Rev. E 74, 056603 (2006). 
T. Kohmoto, H. Tanaka, S. Furue, K. Nakayama, M. Ku- 
nitomo, and Y. Fukuda, Pliys. Rev. A 72, 025802 (2005). 
Md.A.I. Talukder, Y. Amagishi, and M. Tomita, Phys. 
Rev. Lett. 86 3546 (2001). 

A.I. Talukder and M. Tomita, Phys. Rev. A 72, 
051802(R) (2005). 

K.E. Oughstun and N.A. Cartwright, J. Opt. A: Pure 
Appl. Opt. 4, S125 (2002). 

J. Peatross, S.A. Glasgow, and M. Ware, Phys. Rev. Lett. 
84, 2370 (2000). 

J. Peatross, M. Ware, and S.A. Glasgow, , J. Opt. Soc. 
of Am. A 18, 1719-1725 (2001). 

K.E. Oughstun and N.A. Cartwright, J. Mod. Opt. 52, 
1089 (2005). 

Lipsa Nanda, Aakash Basu, and S.A. Ramakrishna, 
Phys. Rev. E 74, 036601 (2006). 

L. Nanda, H. Wanare, and S.A. Ramakrishna, Phys. Rev. 
A 79, 041806(R) (2009). 

A. Kuzmich, A. Dogariu, L.J. Wang, P.W. Milonni, and 

R.Y. Chiao, Phys. Rev. Lett. 86 3925 (2001). 

S.-Y. Zhu, L.-G.Wang, N.-H. Liu, and Y.-P. Yang, Eur. 



Phys. J. D 36, 129 (2005). 
[15] L.J. Wang, A. Kuzmich, and A. Dogariu, Nature (Lon- 
don) 406, 277 (2000). 
[16] S. Chu and S. Wong, Phys. Rev. Lett. 48, 738 (1982). 
[17] M.O. Scully and M.S. Zubairy, Quantum Optics (Cam- 
bridge University Press, Cambridge, U.K., 1997). 
E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou and 
CM. Soukoulis, Nature 423, 604 (2003). 
V. Kuzmiak and A.A. Maradudin, Phys. Rev. B 55, 7427 
(1997). 

M.E. Ta§gm, O.E. Miistecaplioglu and M.O. Oktel, Pys. 
Rev. A 75, 063627 (2007). 

J.D. Jackson, Classical Electrodynamics (Wiley, New 



[18: 

[19 
[20 
[21 
[22; 



[23' 
[24 

[25; 



[26 
[27 



York, 3"'' Ed. 1998) Sees 7.8 and 7.9. 

M. Tanaka, M. Fujiwara, and H. Ikegami, Phys. Rev. A 

34, 4851 (1986). 

R.L. Smith, Am. J. Phys. 38, 978 (1970). 

J.D. Jackson, Classical Electrodynamics (Wiley, New 

York, 3"'^ Ed. 1998) p.s 331,337. 

M.J. Corinthios, lEE Proc.-Vis. Image Signal Process. 
150, 69 (2003); M.J. Corinthios, lEE Proc.-Vis. Image 
Signal Process. 152, 97 (2005). 

We perform the a ^ iz transformation in the Ref. [25l ]. 
Definition of ^(z) = limr^oo /J*^ e~'"* = 
limi— 2sin{zT)/z results no pole. 



[28] |http: / /www.nature.com / news / neutrino-experiment-replicates- faster- 



